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Abstract 

In the quantization scheme which weakens the hermiticity of a Hamiltonian to its 
mere VT invariance the superposition V[x) = + Ze^/x of the harmonic and 
Coulomb potentials is defined at the purely imaginary effective charges [Ze^ = if) 
and regularized by a purely imaginary shift of x. This model is quasi-exactly solvable: 
We show that at each excited, (A^ + 1)— st harmonic-oscillator energy E = 2N + 3 
there exists not only the well known harmonic oscillator bound state (at the vanish- 
ing charge / = 0) but also a normalizable (A^ + 1)— plet of the further elementary 
Sturmian eigenstates ipin}^^) at eigencharges / = /{„} > 0, n = 0, 1, . . . , A^. Beyond 
the smallest multiplicities A^ we recommend their perturbative construction. 



PACS 03.65.Ge, 03.65.Fd 



1 Introduction 

Schrodinger equation for all the asymptotically harmonic oscillators in one dimension, 



need not necessarily be kept defined just on the real axis of coordinates. Indeed, its 
available physical asymptotic solutions 



are normalizable not only on the real intervals x > ^ 1 and x < —Xa but also in 
their complex vicinity defined by the respective formulae 

G (xa, cxd), arg x G (— 7r/4, 7r/4) and arg x G (— 57r/4, — 37r/4). (3) 

One may connect these two complex wedges by a contour of integration which is 
arbitrarily deformed within the domain of analyticity of the potential in question. 
In particular, one may choose the purely harmonic V{x) = x^ and move the real 
axis of coordinates x to a parallel line at a distance a. The new, shifted potential 
V{x) = x^ + 2aix — a? loses its hermiticity and preserves only a certain symmetry 
with respect to a simultaneous change of the parity P (x — > — x) and of the time 
ordering T (which means the mere complex conjugation i ^ ~i'm time-independent 
cases), without any change in the discrete spectrum itself. 

Similar paradoxes seem to have inspired a deeper analysis of the various VT 
symmetric phenomenological models. Bessis and Zinn- Justin related some of 
them to the so called Lee- Yang zeros in field theory 0] and Bender with coauthors 
paid attention to their possible role in the parity breaking 0, phase transitions |^ 
and quantum electrodynamics 

VT symmetry does not prove less exciting on a purely methodical level. Its use 
may range from the fundamental problems of the ambiguities of quantization |p 
and of the exact solvability and Darboux transformations up to the questions 
of convergence of perturbation expansions [|]. In this context, there also appeared 




(1) 



%Ij{x) ~ exp 



2 



ax + 61n (x) + 0(l/x) 



h={E + a'-l)/2 (2) 
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an amazing discovery [^] of the so called quasi-exact (which means incomplete [|T0|) 
solvability of quartic oscillators. 

In contrast to its unsolvable three-dimensional counterpart the latter quartic 
model did not contain the Coulombic component e^/x [|1^. In all the one-dimensional 
Hermitean models this term is routinely being omitted due to its strongly singular 
character in the origin. In the PT— symmetric non-hermitean setting, nevertheless, 
a purely imaginary shift could regularize this singular term in principle. This was 
our main inspiration. 

Even after the screening x —* x—ic the persistent imaginary part of forces V{x) ~ 
l/{x — ic) would still violate unitarity and cause the predominance of unstable, 
resonant bound states. At this point we may recall the above-mentioned strategy 
which tries to compensate the instabilities by the constraint of VT symmetry. This 
re-defines the electric charge. With real / in its effective PT— symmetric value 
Ze'^ = if the nontrivial interaction model need not even contain the cubic and 
quartic terms. Its Schrodinger equation with the three real parameters a, f and c 
reads 

\ 2 f 
— — + X + 2iax + I 

dx-^ X — ic 

and gives, presumably, a real and discrete spectrum of energies E. 

In what follows, we shall re-examine eq. @) from the point of view of its possible 
quasi-exact solvability. We shall be able to show that its elementary (A^ + 1)— plets 
of Sturmian eigenstates exist at any = 0, 1, . . . (Section 2) and that their explicit 
construction (mediated by the vanishing of an (A^ + 1) x (A^ + 1)— dimensional secular 
determinant, see below) may be significantly facilitated via perturbative techniques 
(Section 3). A summary of our results will be outlined in Section 4. 



ip{x) = Eij{x), ij{±oo) = 0, (4) 



2 Quasi-exact solutions 
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2.1 Taylor series and its termination 

Equations (|I]) and (|]) clarify the structure of all the normalizable solutions iplx) of 
eq. (^) near x = ±00. In the light of the identity l/{x — ic) = {x + ic) / {x"^ + c^) our 
model remains smooth and regular at all the finite coordinates and everywhere off the 
complex pole of V{x) at x = ic. In the vicinity of this complex point it is convenient 
to demand that our wave functions vanish, iIj{x) ~ x — ic. Having in mind the 
possible deformations of the contour of integration, such a requirement is immediately 
inspired by the universal and widely accepted regularization of Schrodinger equations 
near their strong singularities ||12[. On this basis let us now try to solve our present 
Schrodinger bound state problem by means of the following elementary harmonic- 
oscillator-like ansatz 



ip^x) = (c + ix) exp 



iax 

2 



N 



(p{x), >^{x) = ^ hn{ixY . (5) 



n=0 



Such a terminating Taylor-series assumption fixes the energy (cf. eq. (^), 

E = 2N + + (6) 

Its insertion in our differential Schrodinger bound state problem leads to the 
equivalent + 1 recurrence relations 

AnK-l + Bnhn + C„/l„+i + Dnhn+2 = 0, n = 0, 1, . . . , A^. 

Keeping in mind that h^i = h^+i = /iiv+2 = we easily derive the values of the 
coefficients, 

Dn = c{n + l){n + 2), C„ = {n + l){n + 2 - lac), 
Bn = -2a{n + 1) - 2c{N +l-n) - f, A„ = -2(A^ +l-n). 
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They may be arranged in a square matrix with four diagonals, 



( Bo Co Do 
A, B, Ci D, 
A2 ••• ■•■ 



D 



iV-2 



Cn-1 



\ 



B 



N ) 



In notation with the row vectors h — {ho, hi, ... , h^) our recurrences may be then 
re-interpreted as a non-hermitean matrix problem 

QP' = 0. 

Normalization hN — i implies its immediate compact and unique solution 



h 



N-k-l 



det 



V 



BN-k C]Sf-k D]s[-k 

-DiV-2 

-Bat-i Cat-i 
An Bn 



(7) 



for all the relevant indices A; = 0, 1, ... , N—1. At the remaining, "redundant" value of 
k — N the left-hand-side symbol vanishes identically since, due to our assumptions, 
h-i = 0. With the linear charge- dependence in Bn = Bn{f) = -Bn(O) — / and in 
Q = Qif) = Q(0) ~ this is equivalent to the current secular equation 



det [Q(0) - //] = 0. 



(8) 



As a polynomial constraint of the (A^^-l-l)— st degree this determines an (A^ + 1)— plet 
of eigencouplings / = f{k}- They are complex in general. In the spirit of our 
introductory remarks it only remains for us to demonstrate that all these roots are 
real (i.e., not breaking the overall T'T— symmetry of our model), in a certain sub- 
domain of parameters at least. 
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2.2 Full VT symmetry in the large-screening domain 



For convenience let us shift the variable / = X + [N + 2)(c — a) and look at Table 1. 
It lists the explicit secular equations at the first few lowest integers A^. From the 
Table we may infer that all the zeros X{n} and/or of eq. are just functions of 
the single parameter = c + a. This is not surprising. The change of the value of a is 
just a shift of the integration path within the domain of analyticity of our potential. 
We may demonstrate the explicit form of this type of invariance of our differential eq. 
(^) algebraically: The equation remains the same after we compensate the change 
of the coordinate x x + i5 by the simultaneous shift a — > a + 5 and c — > c — 5 of 
our pair of parameters. The values of / and d remain unchanged. In our subsequent 
considerations we shall put a = without any loss of generality, therefore. 

Assuming that the strong Coulomb singularity lies off the real line, c 7^ 0, we 
may introduce A = 1/c and re-scale Y = X/c. Our secular equation (||) may be 
re-phrased as corresponding to the asymmetric linear algebraic problem 



^ -r-iv 

-2iVA 



2A 2 
-Y-N + 2 6A 
(2-2iV)A ■•• 



-4A 



V 



N{N~1) 
-Y + N-2 A^(A^+1)A 
-2A -Y + N J 



( 



\ 



ho 
hi 
h2 



hN-i 
\ hN J 



(9) 



In the large-displacement limit A ^ we get an exactly solvable case which deter- 
mines the real eigencharges f{n}- This may be proved easily since in this limit our 
"Hamiltonian" XQ = H — Y I becomes an upper triangular matrix. This leads to 
the following closed formula 

(y|o}, . . . , V}) ^ Wl' • • • ' ^{a!}) = (iV, at - 2, . . . , -at + 2, -N) 

which determines the real and approximately equidistant spectrum of eigencharges 
/ = /{„}(c) = 2{n + l)c + 0{l/c) with n = 0,l,...,N. 
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3 Sturmian bound states at the general N 

Once we fix the integer and recall the nonlinear algebraic definition of our 
(A^+1) eigencouplings / = /{„} we immediately imagine that the explicit construction 
of our Sturmians becomes more and more complicated for the larger dimensions N. 
Even if we extend, accordingly, the list of our secular polynomials in Table 1 we must 
necessarily resort to the purely numerical methods when trying to determine their 
exact roots /. The task is easy for the first few integers only. 

At the large A^, our numerical algorithms may still start from the above equidis- 
tant asymptotic estimates and make use of the expected smooth change of the roots 
/ with the decrease of the parameter \d\ < oo. Indeed, the boundary of the domain 
d > dcriUcai{N) where all our roots / stay real is only slowly growing with A^ since 

dcritical{l) = 2, dcriticalC^) ~ 2.9865, dcritical{^) ~ 3.765 CtC. 

All this indicates that in practice a perturbative evaluation of the eigencharges 
could prove almost as efficient as their direct numerical determination, especially 
beyond the smallest A^. 

3.1 Perturbation expansions with A < XcriUcai 
The inspection of the (extended) Table 1 reveals that 



This means that all the odd perturbation corrections vanish identically, yP^+i] = Q. 
In the similar spirit we may also write 



In order to appreciate the possible merits of such an approach let us return to our 
asymmetric eq. (|^) and notice that at A = it decays into a direct sum of the two 
linear equations. Each of them couples only h^s with the same parity of the subscript 
n. This is a pleasant simplification. For example, the five-dimensional unperturbed 
eigenvalue problem at A^ = 4 decays into the separate two- and three-dimensional 



F(A) = FM + A2r[2] + ;,4y[4]^___ _ 



(10) 



K = /I(a) = /Im + a/Iw + a2/I[2] + ... . 



(11) 
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sub-equations. We may omit the redundant superscripts l*^' and display the latter 
subset for illustration, 

/ 



V 



h2 



0. 



(12) 



-F - 4 2 
-Y 12 
-F + 4 

Relaxing our above normalization convention and working in the integer arithmetics 
(i.e., in full precision, without any round-off errors) its first solution hs^a} = (1,0,0) 
is found for the eigencharge Y = Y^a} = —4 while h{b} = (1, 2, 0) is obtained at the 
vanishing Y^i,} = and h^c} = (3, 12,4) results from Y^^} = +4. 

Returning to our matrix form Q = Oof the differential Schrodinger eq. (Q) with 
"Hamiltonian" \Q = H — Y I and decomposition H = H^'^^ + \ H^^^ at a general N we 
may now solve it by means of the textbook Rayleigh-Schrodinger perturbation theory 



13|. In the present implementation of this recipe the k—th unknown perturbation 



corrections will be determined by the relation 



y[2] (Jl[k-2] 



y w (Jim 



0. 

We have to solve it step-by-step, at all the subsequent perturbation orders C(A^) 
numbered by the integer k = 1,2,.... For illustration let us pick up = 2 and 
y[o] — 2. In the first-order approximation with A; = 1 we normalize h^^ = and drop 
all the superscripts t^' from yt^l and h^^^. The perturbed equation (^ then acquires 
its 0{X) first-order form 

/ 



/ 



V 



-4 





+2 
-2 




V y 



V 



-Y 


2 





-4 


-Y 


6 





-2 


-Y 







0. 



This implies that Y = Y^^'^ = (as expected) and h = h^^^ = (0,4,0). By the way, 
the latter result nicely illustrates a general feature: We might replace eq. ( [Til) by a 
component-by-component expansion in the powers of squares in a way paralleling 
the expansion of charges Y above. This reflects a symmetry of our model with respect 
to the change of sign of the parameter of screening d. 
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3.2 The role of asymmetry 



We have established that our recurrences (P) provide a new, not entirely standard 
framework for application of perturbation theory. The manifestly non-hermitean 
two-diagonal structure of the underlying exactly solvable unperturbed Hamiltonians 
requires a modification of the formalism itself. The necessity of re- considering 
the standard concepts (say, of the model-space projector) seems to deserve a separate 
comment. The point is that due to the manifest difference between our original and 
transposed pseudo-Hamiltonians we must consider not only the "direct" zero-order 
eigenvalue problem [Q{0) — f I] (h^^^^ = (or rather equation 

T 







at a fixed subscript a = 0, 1, . . . A^) but also its transposed, non-equivalent pendant 



(^ti>) 



T 



0. 



It is worth noticing that in a way paralleling eq. (0) above, all the left eigenvectors 
g = g{X) may be defined by closed formulae, 

/ 



dk+l 



{N-k-l)\ 
2^+1 (A; + 1)! 



det 



Bo Co 



Da 



A 



k-l 



V 



Bk-1 Ck-i 
Ak Bk 



Their knowledge would simplify the large-order algorithm. Still, due to the un- 
physical, auxiliary character of all the left eigenvectors, it is shorter to generate 
and use them just in the zero order. Dropping their superscript t''' as redundant, 
we may return to our illustration ([1^) and complement it by the left eigenvectors 
g{a} = (4,-2,3), g{h} = (0,1,-3) and g^c] = (0,0,1). Incidentally, the overlaps 
G = {g, h) form a diagonal matrix, Gu = G33 = 4, G22 = 2. This trivializes 
F = G^^ needed in the model-space projectors P^a} = h'^a}^a,ag{a}- They are, 
counter-intuitively, non-diagonal but compatible with the usual completeness rela- 
tion / = J2a,(3 hJ^}Fa,l3g{p}- 
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4 Concluding remarks 



4.1 Sturmians and the algebra s/(2). 



Our wave function (^) becomes purely real after the change of variables x = —iy 
which rotates the axes by 7r/2. Of course, the exponential growth of our asymptotics 
(§) in the new artificial variable y has no meaning at all and only the original coor- 
dinate X remains physical. Still, the use of y simplifies our ansatz (^). With a = it 
leads to a new differential Schrodinger equation H0{y) = f 0{y) for real polynomials 
0{y) = ip{—iy) themselves. In this notation our third Hamiltonian-like operator H 
is the quadratic function of generators of the complex s/(2) Lie algebra, 

J'^J- + - 2cJ^° + {N + 2) J- - 2{N + c) 



H 



with 



dy' 



J- J' 



d 



J- 



2J^ = 2{y + c)—-2N, 

{y + cf^-2{y + c)N. 
dy 



This fits the general scheme [|14|, parallels closely the similar quartic-oscillator result 
of ref. [0 (with the same algebra but different H) and re-confirms our above conclu- 
sion that in spite of its non-hermiticity, our PT— symmetric Schrodinger eq. @) is 
quasi-exactly solvable. 

In this context, we would like to emphasize that there exists a nice parallel be- 
tween the Hermitean and PT— symmetric quasi-exactly solvable models. Picking 
up the characteristic examples and adding our present results, we may now sum- 
marize and list the following four different possibilities and types of the quasi-exact 
constructions. 

• The characteristic polynomial example V{x) = a x"^ + (3 x'^ + x^ with a con- 
strained variability, say, of the coupling a leads to the existence of some 
elementary bound states at a finite multiplet of the binding energies in the 
hermitean case [TH . 
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• With a constrained variability of the energy E the non-polynomial hermitean 
potentials exemplified by V{r) = + a/(l + (3 x"^) lead to the elementary 
Sturmian solutions at a finite A^— plet of couplings a |]16| . 

• In the "PT— symmetric case, the complex quartic polynomial potential of ref. 
1^, V{x) = a X + (3 x"^ + 'J — x'^ with a constrained variability of a, exhibits 
the quasi-exact solvability of the sextic-oscillator type. 

• V{x) = x'^ + ia x + i(3 (a: + ?7)/ (0:^ + 7^) of the present model (^) is the "missing" 
PT— symmetric partner to the quasi-exact solvability of the non-polynomial 
Sturmian type. 



4.2 Complex charges 

We have seen that the one-dimensional Schrodinger equation becomes quasi-exactly 
solvable for Coulomb plus harmonic superpositions of potentials provided only that 
we regularize these forces in the "PT— symmetric manner. Multiplets of the Sturmian 
eigenstates acquire then an elementary polynomial form at certain purely imaginary 
couplings Ze^ = i /{„} at = 0, 1, . . . , and any = 0, 1, . . .. 

In a tentative physical support of these complex electric charges we might recall 
a few of their VT symmetric predecessors. A close connection exists with the Bessis' 
cubic force possessing a purely imaginary coupling g = i f. Its rigorous mathematical 
analysis has already been delivered, many years ago, by Calicetti et al . One may 
also mention an even closer parallel with the Bender's and Milton's electrodynamics 
which replaces, for several independent reasons |p, the charge e itself (i.e., not its 
present square e^) by a purely imaginary quantity. 

In a broader methodical context, our new phenomenological model (D extends 
the family of potentials which are comparatively easily described in the language of 
analytic continuations |TB[. Its distinctive feature, in this setting, is the presence of 
a complex pole. Its most important formal merit is its quasi-exact solvability shared 
with the quartic model of ref. H]. In the latter comparison, one could emphasize the 
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"more natural" asymptotic behaviour of our present wave functions: The real axis of 
coordinates still lies within the wedges of the admissible analytic continuations. 

In the context of perturbative considerations it is amusing to notice that a "nice" 
(which means real, discrete and bounded) character of spectrum of our present VT 
symmetric model may be viewed as a result of perturbation of any one of its two 
exactly solvable halves. 
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Table 1. The first six secular equations (H) for eigencharges f +{N +2){a—c) = X, 
with a variable parameter c + a = d and abbreviation d"^ — N — 3 = h. 



N 


detg(/) = 







-X = 


1 




X^-h = 


2 




-X^ + AhX + 8d = 


3 




10 hX^ - 48 + 9 /i^ - 36 = 


4 


-X^ + 20 hX^ 


+ 168 dX^ - 32 (2 - 9) X - 384 hd = 


5 


X^ - 35 hX^ - 


448 dX^ + (259 - 1296) X'^+ 






+3520 hdX - 225 + 10000 h + 51200 = 
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